Detection of Large Scale Structure in a B < 17 m Galaxy Redshift Survey. 



On 
On 

A. Ratcliffe*, T. Shanks*, A. Broadbent*, Q.A. Parkerj, 
^ ! F.G. Watsonf, A.P. OatesJ, R. Fong* & C.A. Collins*; 



*Physics Dept., Univ. of Durham, South Road, Durham, DH1 3LE, U.K. 
| Anglo- Australian Observatory, Coonabarabran, NSW 2357, Australia. 



> 

. \Royal Greenwich Observatory, Madingley Road, Cambridge, CBS OHA, U.K. 

o 

(N 

^Astrophysics Group, Liverpool- J ohn-Moores Univ., Liverpool, L3 3AF, U.K. 
We report on results from the Durham/UKST Galaxy Redshift Survey where 

o 

H 

we have found large scale "cellular" features in the galaxy distribution. These 

cd ' 

have spatial 2-point correlation function power significantly in excess of the 

• i-H 

5_i ■ predictions of the standard cold dark matter cosmological model 1 , supporting 

the previous observational results from the APM survey 2 ' 3 . At smaller scales, 
the 1-D pairwise galaxy velocity dispersion is measured to be 387lg2 kms -1 
which is also inconsistent with the prediction of the standard cold dark matter 
model 1 . Finally, the survey has produced the most significant detection yet of 
large scale redshift space distortions due to dynamical infall of galaxies 4 . An 
estimate of fi° 6 /b = 0.55 ± 0.12 is obtained which is consistent either with a 
low density Universe or a critical density Universe where galaxies are biased 
tracers of the mass. 

The principal aims of the Durham/UKST Galaxy Redshift Survey are to investigate 



the structure and dynamics of the Universe on scales l-100fo _1 Mpc (where h is Hubble's 
constant in units of 100 kms -1 Mpc -1 ). The survey was constructed by measuring the 
redshifts for 1 in 3 of the galaxies from the Edinburgh/Durham Southern Galaxy Catalogue 5 
(EDSGC) to bj ~ 17 m . The resulting survey contains ~ 2500 redshifts, covers a ~ 20° x 75° 
contiguous area of the sky at the South Galactic Pole and probes to a depth of > 300/t -1 Mpc 
with a median depth of ~ 150/i _1 Mpc. The total volume of space surveyed is ~ 4 x 
10 6 /i _3 Mpc 3 , approximately six times that of the updated Southern Sky Redshift Survey 6 , 
due to our ~ 2.5 m fainter limit. 

Figure 1 shows the distribution of the galaxies in right ascension and velocity for each 
of four, 5°, declination slices. The survey gives the striking impression that the galaxy 
distribution is "cellular" or "bubble-like" on 50-100ft, _1 Mpc scales. The most noticeable 
structure is the low density region lying between and 90fo _1 Mpc surrounded by long 
"walls" of galaxies. This structure is present in the slices at —30°, —35° and —40° and 
was referred to as the Sculptor Void in earlier, shallower surveys 7 ' 8 . In the most northerly 
slice at —25° there is evidence that we are seeing the top of this structure and that it 
is indeed a "cell". Figure 2 shows the histogram of galaxy number with distance which 
has several large peaks, also indicating the presence of large scale structure. However, it 
should be noted that at least one of these peaks (eg. at ~ 90/i~ 1 Mpc) does not follow the 
128/i _1 Mpc periodic pattern previously claimed 9 along the North-South Galactic Pole axis 
which intersects our survey at ~ h 54 m , -27.5°, crossing the nearer and further Sculptor 
superclusters. Figure 1 clearly reveals a galaxy distribution that is more complex than any 
simple periodic pattern. 

To test the significance of these large scale features, we have made a statistical investi- 



gation of the galaxy distribution in the Durham/UKST survey by calculating the redshift 
space 2-point correlation function 10 , £(s). This is a measure of the excess probability of 
finding two galaxies in given volume elements at separations s over what is expected for a 
homogeneous galaxy distribution. In figure 3a we compare the Durham/UKST £(s) with 
that from other optical redshift surveys 3,11,12 . On small scales (< 10/i _1 Mpc) we see that 
there is good agreement between all the surveys. On large scales (> 107t _1 Mpc) we see 
that our new estimate is consistent with the detection of large scale power from the APM- 
Stromlo survey but less consistent with the previous Durham surveys. This inconsistency 
is thought to be partly statistical and partly due to the previous use of a simple unweighted 
estimator which, although having minimum errors, was systematically biased against the 
detection of large scale power (Ratcliffe et al. , in preparation) . Our new estimate is also 
consistent with that found from the QDOT infrared redshift survey 13 . Figure 3b shows the 
Durham/UKST £(s) with the results from mock catalogues drawn from two sets of cold 
dark matter 1 (CDM) N-body simulations 14,15,16 ; standard CDM with Qh = 0.5, b = 1.6 
(SCDM) and CDM with Qh = 0.2, 6 = 1 and a cosmological constant (A ^ 0) to ensure 
a spatially flat cosmology (LCDM), where f2 is the mean mass density of the Universe 
and b is the linear bias parameter relating the galaxy and underlying matter distributions, 
(Ap/p) g = b(Ap/p) m . The mock catalogues were constructed with the same angular and 
radial selection functions as the Durham/UKST survey. On small scales (< 10fo _1 Mpc) 
both CDM models agree well with the data. On large scales (> 10/i _1 Mpc) SCDM shows 
no significant large scale power whereas LCDM shows significant power out to ~ 30ft, _1 Mpc. 
The data therefore shows significant excess power at ~ 15-30/i _1 Mpc over SCDM (> 3a). 
LCDM is more consistent with the data at these scales, although even this model produces 



too little power at the 2a level. This rejection of SCDM in our new survey is consistent 
with the findings from the APM correlation functions 2 ' 3 . 

Inferring galaxy distances from redshifts results in the distortion of the clustering pattern 
because of the non-Hubble component of galaxy velocities. This anisotropy along the line 
of sight can be used to estimate some important cosmological parameters describing the 
dynamics of the Universe 4 ' 10 . Figure 4a shows a contour plot of the 2-point correlation 
function, £ v (a, 7r), where a and ir are the perpendicular and parallel separations to the line 
of sight, respectively. On small, non-linear scales the contours are elongated parallel to the 
line of sight caused by the rms virial velocities of galaxies. Using standard techniques 10 to 
model this smearing of £ by the peculiar velocities, we find < v 2 >2= 3871 j£j kms^ 1 at 
perpendicular separations a < 2/i~ 1 Mpc. Assuming that simple biasing models apply 17 , 
this is inconsistent with the SCDM prediction 1 of ~ 1000 kms -1 but nearer the LCDM 
prediction 1 of ~ 600 kms^ 1 . On large, linear scales the contours are compressed parallel 
to the line of sight caused by the infall of galaxies into overdense regions. We use two 
independent methods to estimate Q os /b as shown in figure 4b. Firstly, using the spherical 
harmonic moments 18 (£) of £ v (cr, n) to measure the shape of the contours of constant £, we 
find ri°- 6 /fr = 0.55 ± 0.12 in the region 9-23/i _1 Mpc. Secondly, using the enhancement in 
redshift space 4 of the volume integral of £ (J3) over the real space one, we find Q os /b = 
0.48 ± 0.18, quoting the result at 10/i _1 Mpc, which is representative. Our measured values 
of f2°' 6 /6 can be compared with other optical values of Q°' G /b estimated from redshift space 
distortions, namely 0.48 ± 0.12 (Loveday et al. , manuscript submitted), 0.77 ± 0.16 19 and 
0.5 ± 0.25 20 . All of these values agree within 2a of our results. Thus taking two fiducial 
values for b, b = 1 implies 0, = 0.37 ± 0.15 and 6 = 2 implies Q = 1.17 ± 0.31. This result 



tends to argue against an unbiased critical density, fi = 1, Universe, but favours a biased, 
S7 = 1, Universe or an unbiased, low J7, Universe. 

Thus our Sl 6 /& result agrees with the SCDM and LCDM models but, since they both 
predict S7 6 /6 ~ 0.4-0.6, we cannot discriminate between them. However, we have shown 
that SCDM underpredicts our 2-point correlation function results at large scales and over- 
predicts the 1-D pairwise velocity dispersion at small scales. Therefore, overall, our results 
argue for a model with a density perturbation spectrum more skewed towards large scales, 
such as LCDM. 
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Figure 1. 

Redshift cone-plots of the galaxy distribution in the four contiguous 5° declination strips 
from the Durham/UKST redshift survey, (a), (b), (c) and (d), respectively. Large "cellular" 
structures in the galaxy distribution on scales of 50-100/i~ 1 Mpc are clearly seen. The 
redshifts of these galaxies were obtained over the three year period 1991-1994 using the 1.2m 
UK Schmidt Telescope at Siding Spring, Australia, in conjunction with the fibre-coupled 
spectroscopy system FLAIR 21 . The completeness is > 75% to the nominal magnitude limit 
of bj = 17.0 m . Our magnitude limit varies slightly (0.1 m -0.2 m ) from field to field but this 
is exactly taken into account in all our statistical analysis. A comparison of 150 galaxies 
with published velocities 7,8 ' 22 ' 23 indicates that our measurements have a negligible offset 
and are accurate to ±150 kms^ 1 (Ratcliffe et al. , in preparation). 

Figure 2. 

The histogram of galaxy number with distance, n(r), in the Durham/UKST survey. The 
dashed curve shows how a random and homogeneous distribution would appear given our 
angular /radial selection functions, sampling rate and galaxy luminosity function from the 
Durham/UKST survey (Ratcliffe et al. , in preparation). The arrows indicate where the 
128/i _1 Mpc "spikes" in the galaxy distribution 9 should appear in this region of the sky. 
In this larger angle survey the galaxy distribution clearly appears more complex than any 
simple periodic pattern. There are two strong peaks at ~ 90 and ~ 170/i _1 Mpc signifying 
"walls" in the galaxy distribution. There is possible evidence for a third such feature at 
~ 270/i- 1 Mpc. 



Figure 3. 

In (a) we compare the Durham/UKST 2-point correlation function, £(s), with that from 
other optical redshift surveys 3,11,12 , £(s), was calculated using the weighted estimator 24 ' 25 
which gave no systematic offset and the minimum variance in the CDM mock catalogues. 
On this plot the error bars for the Durham/UKST survey were calculated using the la 
variance from the LCDM mock catalogues. In (b) we compare with a SCDM model and 
a LCDM model. The shaded areas are the la confidence regions on an individual mock 
catalogue for the above two models of CDM using the aforementioned optimal estimator 
and weighting combination so that a direct comparison between the data and models can be 
made. We have also estimated the integral constraint to see if it would cause any systematic 
offset in £(s) on large scales and find this to be a negligible effect. 

Figure 4. 

In (a) we show the contours of constant 2-point correlation function where the effects of 
small and large scale redshift space distortions are clearly seen. Solid lines denote £ > 1 
with A£ = 1, short dashed lines denote < £ < 1 with A£ = 0.1 and long dashed lines 
denote £ < with A£ = 0.1. For reference, the lines with £ = 1 and £ = are in thick 
bold and the regularly spaced thin bold lines show an isotropic correlation function for 
comparison. A representative error bar is also plotted. In this plot £„(<7, n) was calculated 
using a single pair weighting of galaxies 10 because this weighting produced less noise in 
£u(<7, 7r). Our consistent use of this different weighting scheme has no effect on the either 
method of estimating Q 0M /b (Ratcliffe, et al. in preparation). In (b) we show our estimated 
values of Q os /b for two methods 18 ' 4 as a function of separation r. The solid line shows 



the maximum likelihood fit for Q°- 6 /b to the first method of estimation (£) in the region 
9-23/i _1 Mpc and the dashed lines denote the la confidence intervals of this fit. The error 
bars on these independent points come from the LCDM mock catalogues. The error bar of 
the estimated point at ~ 10/i _1 Mpc from the second method (J3) is representative of the 
errors from this method and was calculated from the variance between the 4 quadrants of 
the survey. Only one error bar is shown because of the non-independence of the points in 
this method. Both of these sets of error bars include the contribution from cosmic variance. 
We have estimated £ in real space using two methods involving the inversion of the projected 
correlation function 26 (Ratcliffe et al. , in preparation) and find consistent results. 
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